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SLAB  MODEL  ANALYSIS  OF  MAGNETIC  AND  COLLISIONAL  VISCOSITY 
EFFECTS  ON  THE  FIRST  GENERATION  OF  NUCLEAR  STRUCTURE 


I.  INTRODUCTION 

Our  goal  in  this  paper  is  to  postulate  a  simplified  model  of  the 
structure  formation  process  in  the  late  time  nuclear  environment.  To  the 
extent  that  this  model  is  valid,  the  specification  of  the  freezing  scale  in 
the  nuclear  case  is  reasonably  simple  both  conceptually  and 
computationally.  As  in  the  case  with  any  theory,  its  success  depends 
critically  on  its  assumptions,  and  ve  ask  the  reader  to  examine  our 
hypotheses  with  a  critical  eye. 

In  effect,  we  are  going  to  argue  that  we  can  use  one-dimensional 
linear  stability  analysis  to  predict  the  freezing  scale  in  the  late-time 
mul t i d i mens i onal  nonlinear  evolution  of  a  nuclear  plume.  In  this  respect 

our  view  is  quite  similar  to  that  of  the  JAYCOR  freezing  model  [Glassman 

and  Sperling,  1987].  The  present  arguments  differ  from  those  used  by 
JAYCOR  in  two  key  respects:  (1)  the  success  of  the  model  relies  on  the 
shielding  properties  of  multidimensional  high  M  clouds,  rather  than  a 
hoped-for  proportionality  between  one-and  two-dimensional  stability 
analysis  scale  sizes;  and  (2)  we  find  that  the  long  wavelength  limit  is 

often  not  appropriate  for  the  range  of  nuclear  parameters  we  expect.  Only 

this  last  point  has  any  effect  on  the  model  itself,  and  it  is  fair  to  say 
that  the  only  structural  difference  between  the  model  we  propose  here  and 
that  of  JAYCOR  is  our  use  of  routines  which  solve  the  full  stability 
problem  for  the  actual  expected  density  gradients,  rather  than  for  the 
asymptotic  long  wavelength  limit.  Thus  changing  a  piece  of  software  which 
implements  the  present  JAYCOR  model  to  the  model  ve  propose  here  would 
simply  involve  changing  the  routines  which  calculate  growth  rates  to  our 
routines,  or  better  yet,  to  a  table  look-up  based  on  a  one-time  parameter 
search  with  our  routines. 

Miauacripc  approved  October  20.  1988. 
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In  concise  form  our  postulated  model  is  as  follows: 

1.  The  one-dimensional  slab  mode  is  appropriate  for  the  linear  regime 
during  the  first  generation  of  nuclear  structure  (see  Section  II). 

2.  Due  to  the  high  M  ratio  of  a  nuclear  plume,  and  to  the  scale  sizes 

of  the  neutral  wind  relative  to  k~*  (see  below),  the  first  generation  of 

nuclear  structure  is  the  last  (see  Sec.  III). 

3.  If  k  is  the  wavenumber  for  which  the  growth  rate  y(k)  maximizes, 

m 

and  if  the  maximum  in  y  is  sufficiently  narrow  about  k  ,  then  the  nonlinear 

n 

"blobs"  associated  with  the  first  generation  of  nuclear  structure  are 

dominated  by  a  scale  size  k-1  (see  Section  IV). 

m 

4.  It  follows  from  (1)  -  (3)  that  the  freezing  scale  is  given  by  k'1 

In  Section  II  we  make  the  case  for  the  use  of  one-dimensional  slab 

linear  stability  analysis  to  model  the  behavior  of  a  two  dimensional 

nuclear  plume.  In  Section  III  we  make  the  case  for  the  first  generation  of 

nuclear  structure  being  the  last.  In  Section  IV  we  show  numerically 

generated  curves  of  the  growth  rate  y  for  a  range  of  expected  nuclear 

parameters,  and  show  that  when  the  function  y(k)  has  a  sufficiently  narrow 

maximum  at  wavenumber  k  ,  then  the  structure  scale  size  at  the  end  of  the 

m 

first  generation  of  structure  will  be  dominated  by  k^-  .  Combining  this 
with  the  results  of  Section  II,  we  conclude  that  the  freezing  scale  is 
simply  k  In  Section  V  we  show  the  results  of  some  two-dimensional 

nonlinear  numerical  simulations,  which  verify  what  we  have  concluded  in  the 
previous  sections.  tfe  also  address  in  this  section  the  question  of  the 
sensitivity  of  our  conclusions  to  the  size  of  the  seed  perturbations  that 


2 


T 


are  assumed  to  exist  in  all  real-vorld  nuclear  plumes.  In  Section  VI  we 
state  our  conclusions,  and  try  to  place  this  model  in  perspective  by 
enumerating  those  things  which  may  limit  its  applicability.  The  physics 
and  basic  equations  used  here  are  given  in  Zalesak  et  al.  (1988)  along  with 
the  linearized  equations  used  in  the  stability  analysis. 

II.  THE  CASE  FOR  USING  ONE-DIMENSIONAL  GROWTH  RATES 

Consider  a  nuclear  plume  at  the  beginning  of  the  late-time  environment 
(Fig.  1).  We  assume  that  the  overall  cross-field  diameter  of  the  plume  is 
tens  or  hundreds  of  kilometers,  with  edge  gradient  scale  lengths  in  the 
kilometer  range,  certainly  no  shorter  than  100m,  and  that  the  ratio  M  of 
integrated  Pedersen  conductivity  inside  the  plume  to  outside  the  plume  is 
very  large  (100  or  greater).  The  two  primary  drivers  for  structure  are  (1) 
a  component  of  the  neutral  wind  more  or  less  uniform  and  imposed  on  the 
plume  externally,  and  (2)  those  components  of  the  neutral  wind  which  are 
internal  to  the  cloud,  i.e.  a  radial  neutral  wind  or  at  least  one  which  has 
substantial  variation  over  the  diameter  of  the  plume  (see  Fig.  1). 

Ve  know  that  the  external  component  of  the  wind  will  not  be  effective 
in  producing  structure,  due  to  the  large  M  of  the  cloud.  The  cloud  will 
simply  produce  polarization  electric  fields  ("shielding")  which  allow  it  to 
drift  at  the  field-line-veighted  neutral  wind  velocity,  reducing  the  field- 
line-weighted  slip  velocity  and  driving  the  growth  rates  to  near  zero 
[McDonald  et  al.,  1981).  The  plume  will  steepen  on  one  side  as  a  result  of 
this  process  however.  Thus  we  argue  that  it  is  the  internal  winds  which 
have  the  only  real  chance  of  producing  structure  in  the  nuclear  plume, 
since  shielding  cannot  take  place  in  the  internal  wind  geometry,  and 
further,  that  to  the  degree  that  this  internal  wind  component  is  well 
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approximated  by  a  radial  wind,  the  structuring  process  is  veil  described  by 
a  slab  model  analysis  (see  Fig.  2). 

It  remains  to  specify  the  gradient  scale  lengths  and  slip  velocities 
to  be  used  in  the  slab  analysis.  In  Fig.  1  ve  show  the  range  of 
possibilities  that  ve  consider  to  be  most  likely.  In  Fig.  la  ve  shov  a 
reasonably  symmetric  nuclear  plume  located  near  the  center  of  a  radial 
neutral  wind  field  associated  vith  a  burst.  Since  the  rav  vind  speeds  are 
expected  to  be  quite  high  near  this  point,  and  since  the  radial  geometry 
yields  little  or  no  shielding,  ve  expect  the  slip  velocities  to  be  large. 
Further,  since  no  shielding  means  no  steepening  either,  ve  expect  the 
gradient  scale  lengths  to  be  large.  Assuming  that  structuring  does  not 
take  place  at  this  time,  a  little  vhile  later  ve  expect  the  plume  to  appear 
as  in  Fig.  lb.  The  small  external  component  of  the  neutral  vind  has 
translated  the  entire  cloud  avay  from  the  neutral  vind  center,  reducing  the 
rav  neutral  vind  speed  and  also  its  internal  component,  thus  giving  smaller 
slip  velocities.  Hovever  the  translation  has  also  given  rise  to  a 
steepening  on  one  side,  thus  giving  smaller  gradient  scale  lengths.  Fig. 
lc  just  continues  the  arguments  advanced  to  justify  Fig.  lb,  vith  still 
smaller  slip  velocities  and  gradient  scale  lengths.  Thus  ve  conclude  that 
ve  can  span  the  range  of  realistic  possibili ties  by  pairing  small  gradient 
scale  lengths  vith  small  slip  velocities,  and  large  gradient  scale  lengths 
vith  large  slip  velocities.  Note,  hovever,  that  this  raises  the  distinct 
probability  that  the  long  vavelength  limit  vill  not  alvays  be  appropriate 
for  the  slab  stability  analysis,  and  that  finite  gradient  scale  lengths 
vill  have  to  be  part  of  the  analysis. 
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III.  THE  CASE  FOR  THE  FIRST  GENERATION  OF  NUCLEAR  STRUCTURE  BEING  THE  LAST 

The  above  simply  argues  that  a  slab  analysis  is  sufficient  to  predict 

the  first  generation  of  structure  in  a  nuclear  environment.  In  order  to 
fully  argue  the  case  for  using  slab  analysis  to  predict  the  late-time 
freezing  scale,  ve  must  also  make  the  case  for  the  first  generation  of 
nuclear  structure  being  the  last.  This  ve  do  as  follows. 

Consider  the  cases  shown  in  Fig.  1.  Ve  shall  assume  that  the  scale 
lengths  associated  with  the  divergent  part  of  the  neutral  flow  field  are 
approximately  equal  to  those  of  the  large  unstructured  plume,  on  the  order 
of  tens  or  hundreds  of  kilometers,  since  both  evolved  from  the  same  or  a 
similar  burst.  Thus  when  structuring  takes  place,  generating  kilometer  or 
sub-kilometer  scale  plasma  structures,  the  neutral  wind  field  will  now 
appear  to  be  almost  totally  externally  imposed  to  these  new,  much  smaller 
structures.  Thus  the  structuring  event  itself  will  have  given  rise  to  a 
sudden  and  substantial  drop  in  the  internal  wind  component.  Ve  now  argue 
simply  that  this  sudden  drop  has  put  these  new  small-scale  structures  in  a 
regime  where  even  the  smallest  of  dissipation  mechanisms  will  inhibit 
bifurcation.  Thus  structuring  ceases  after  the  first  generation. 

IV.  LINEAR  ANALYSIS  FOR  NUCLEAR  PARAMETERS 

In  this  section  we  attempt  to  span  the  range  of  possible  nuclear 

parameters  and  compute  the  growth  rate  y(k)  for  all  relevant  values  of 

wavenumber  k,  based  on  the  conclusions  of  the  previous  sections.  Ve  shall 

be  looking  for  the  situations  where  r(k)  is  a  strongly  peaked  function  of 

k,  since  in  this  case  ve  can  argue  that  at  the  end  of  the  first  generation 

of  structure,  the  resultant  nonlinear  evolution  v'll  be  dominated  by  scale 

sizes  of  order  k~S  where  k  is  the  wavenumber  for  which  r(k)  maximizes, 
m  m 

* 
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Rather  than  try  to  present  the  detailed  argument  now,  let  us  present  a  few 
calculations  of  linear  growth  rates  for  nuclear  parameters  and  use  these  as 
models  for  our  arguments.  The  geometry  we  shall  work  with  is  shown  in  Fig. 
2.  We  have  a  driver  modeled  as  a  neutral  wind  blowing  through  the  edge  of 
a  plume.  This  is  meant  to  represent  a  small  pie-slice  of  a  radial  wind 
blowing  outward  into  a  transition  region  between  the  plume  interior 
electron  density  n<  and  the  ambient  electron  density  n>. 

We  must  choose  the  following  parameters  when  performing  a  linear  slab 
analysis. 

1)  n^,  the  electron  density  outside  the  plume, 

2)  n<f  the  electron  density  inside  the  plume, 

3)  the  functional  dependence  n(x)  in  the  transition  between  n>  and 

V 

4)  the  neutral  wind  slip  velocity  U, 

5)  The  ion-neutral  collision  frequency  \>.n, 

6)  The  ion  species  (e.g.  0*,  B  +),  and 

a 

7)  the  ion  temperature  T^. 

For  the  calculations  we  present,  we  have  chosen  the  following  as  being 
representative  of  the  nuclear  plume  case: 

1)  n>  *  105, 

2)  n<  *  108, 

4)  the  ion  species  is  0+, 

3)  n(x)  is  a  hyperbolic  tangent  with  scale  length  L, 

5)  the  ion  temperature  Ti  =0.1  ev, 

6)  three  calculations  were  done,  in  accordance  with  the  arguments 
advanced  in  Section  II: 

L)  large  L,  large  U  :  L  =  1  km,  U  =  1  km/sec 

M)  moderate  L,  moderate  U:  L  100m,  U  =  100  m/sec 
S)  small  L,  small  U  :  L  =  20m,  U  =  20  m/sec, 
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7) 


for  each  of  those  three  calculations,  three  \>^n  values  were  used, 

meant  to  span  the  altitude  range  from  low  to  highs  v.  »  10.0, 
-1  ln 
1.0,  0.1  s  ,  and 

8)  the  physics  included  are  ion-neutral  collisions,  ion  inertia, 
magnetic  viscosity,  and  ion  collisional  viscosity.  A  detailed 
description  of  the  physics  and  the  numerical  solutions  of  the 
linearized  equations  is  found  in  Zalesak  et  al. ,  1988. 

Calculations  L,  M,  and  S  are  shown  in  Figs.  3,  4,  and  5  respectively. 

In  each  plot  is  shown  the  growth  rate  y  as  a  function  of  the  wavenumber  k, 

for  each  of  the  three  values  of  w.  . 

in 

Let  us  pick  a  particular  growth  curve,  that  of  calculation  M  (Fig.  4) 
for  -  1.0.  Suppose  we  assume  that  our  initial  condition  consists  of 
the  slab  itself  plus  a  small  amplitude  white  noise  perturbation.  The  curve 

peaks  at  k  »  k^  -  20  km-1,  with  growth  rates  at  k  =  10  km-1  and  k  =  30  km'1 

lower  by  a  factor  of  about  1/6.  All  these  modes  will  grow,  but  at 
different  rates.  If  we  postulate  that  the  initial  normalized  amplitude  of 
each  mode,  Sn/n,  is  ~  10'*,  then  it  will  take  13.8  growth  times  for  a  mode 
to  reach  Sn/n  -  1.  Now  consider  two  modes  whose  growth  rates  differ  by 

1/6.  By  the  time  the  faster  growing  of  the  two  had  reached  Sn/n  =  1.0,  the 

slower  growing  one  would  have  achieved  only  Sn/n  =  0.1.  Thus  we  would 
expect  the  spectrum  at  the  transition  to  nonlinear  (Sn/n  =  1)  to  be 
strongly  peaked  at  km  *  20  k~  .  If  we  further  argue  that  these  fingers  are 
now  two-dimensional  structures  and  high  M  and  hence  virtually  immune  to 
further  bifurcation,  then  we  conclude  that  the  freezing  scale  is  k'1. 

Note,  however,  that  not  all  of  the  curves  shown  display  this  property 
of  being  strongly  peaked  functions  of  k.  In  general  for  larger  L  and  for 
larger  v^n,  the  maximum  is  less  strongly  peaked,  and  indeed  may  be 
nonexistent.  In  these  cases  we  expect  our  model  to  break  down. 
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V. 


TWO-DIMENSIONAL  NONLINEAR  SIMULATIONS 


In  this  section  we  take  one  of  the  initial  conditions  considered  in 
the  previous  sections,  the  case  of  L  =  100  m,  Un  =  100  m/sec,  =  1.0 
s"1,  which  yields  a  reasonably  peaked  profile  y(k),  and  put  our  ideas  to 
the  test.  We  seed  our  initial  profile  with  a  white  noise,  random  phase 
initial  perturbation  in  the  x  direction,  in  effect  seeding  all  k's  with  the 
same  level  of  fluctuation.  We  then  let  the  perturbed  slab  evolve  in  time 
using  a  fully  two  dimensional  nonlinear  numerical  simulation  code,  to  see 
if  a  preferred  freezing  scale  does  arise  out  of  the  noise  at  late  times. 
The  curve  for  \>^n  »  1.0  s~*  in  Fig.  4  shows  that  the  peak  in  the  growth 
rate  occurs  for  k  -  20  km-^1.  Thus  we  expect  a  wavelength  of  about  300  m  to 
dominate  at  late  times.  The  grid  used  in  the  code  was  1200  m  wide  in  the 
x-direction,  with  periodic  boundary  conditions  assumed.  Thus  we  expect 
approximately  four  major  fingers  of  plasma  to  emerge  from  the  slab  at  late 
times. 

We  have  run  the  calculation  at  two  different  levels  of  fluctuation. 
The  fluctuation  was  applied  by  initially  displacing  the  plasma  a  distance 


5y  in  y,  where 

5y  is  a  white 

noise  random 

phase  function 

of  X. 

calculations  A 

and  B  we  took  the 

mean  Sy 

to 

be  10  ^  L  and 

5  x  10 

respectively. 

Obviously  we  are  trying 

to 

ascertain  the 

effect 

fluctuation  level  on  our  results,  since  the  logic  we  have  used  in  section 
IV  would  argue  that  the  late-time  spectrum  would  be  less  dominated  by  A  ~ 
300  ra  for  the  high  fluctuation  case. 

In  Fig.  6  we  show  the  isodensity  contours  for  the  initial  conditions 
for  calculation  A.  Note  that  the  perturbation  level  is  sufficiently  low  as 
to  make  the  perturbations  almost  undetectable  on  the  plot.  In  Fig.  7  we 
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show  calculation  A  after  10.67  sec  of  evolution.  Note  that  there  are 
indeed  3-4  major  fingers  of  plasma,  as  the  present  model  would  predict. 

In  Fig.  8  we  show  the  initial  conditions  for  calculation  B.  Note 
that  the  perturbation  is  now  easily  visible  to  the  eye.  In  Fig.  9  we  show 
calculation  B  after  3.56  sec  of  evolution.  Despite  the  presence  of  high 
frequency  structure,  the  morphology  is  still  dominated  by  approximately 
four  major  fingers  of  plasma.  Thus  even  when  the  perturbation  levels  are 
reasonably  high,  the  model  advanced  here  may  be  acceptable. 

VI .  CONCLUSIONS 

Ve  have  postulated  a  simple  model  of  the  freezing  scale  in  nuclear 
plumes,  whose  success  depends  strongly  on  the  high  M  ratios  of  nuclear 
plumes,  and  on  the  existence  of  some  mechanism  which  causes  slab  growth 
rates  y(k)  to  be  sharply  peaked  functions  of  k.  Although  the  results 
examined  here  rely  on  the  and  terms  in  Braginskii's  stress  tensor, 
this  need  not  be  the  case.  The  final  model  is  similar  to  the  JAYCOR  model 
in  the  sense  that  the  freezing  scale  is  taken  to  be  proportional  to  k'1 
where  k  is  the  wavenumber  for  which  Y(k)  maximizes.  However,  the 
derivation  is  different  and,  ve  think,  rests  on  more  solid  theoretical 
ground.  It  is  also  clearer,  making  the  limitations  of  the  model  more 
apparent.  Using  numerical  simulations,  we  conclude  that  the  model  has  a 
reasonable  chance  of  working  in  the  nuclear  case,  with  the  caveat  that  the 
long-wavelength  limit  may  not  always  be  the  appropriate  model  for  the  slab 
growth  rates. 
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Fig.  1.  Range  of  possibilities  for  the  morphology  of  a  nuclear  plume 
perpendicular  to  the  magnetic  field,  during  the  late-time  regime. 
The  arrows  denote  net  tral  wind  direction.  The  possibilities  are 
(a)  L  Large,  Ugj^  l?rge,  steepening  small;  (b)  L  moderate,  Uslip 
moderate,  steepening  moderate;  and  (c)  L  small,  small, 

steepening  large.  See  text  for  reasoning. 
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(a) 


PLUME 

INTERIOR 


AMBIENT 

PLASMA 


(b) 

Fig-  2.  Approximating  a  radial  vind  component  acting  outward  on  a  curved 
cloud  edge  segment  (a)  as  a  uniform  vind  acting  on  a  slab  (b). 
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GAMMA ( SEC- I ) 


M  =  1000.,  EL  =  1000. M,  U0  =-1000. M/S 


Fig.  3.  Plots  of  y  vs  k  for  the  case  of  n>  =  10^,  M  =  1000,  Un  =  -1000 
m/s,  using  a  hyperbolic  tangent  density  profile  and  L  =  1000  m. 
Curves  A,  B,  and  C  refer  to  \>^  =  0.1,  1.0,  ..nd  10.0  s-* 

respectively.  (Calculation  L) 
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GAMMA (SEC-1 ) 


”  1  '  T 


M  =  1000.,  EL  =  100. M,  U0  =  -100. M/S 


Fig.  4.  As  in  Fig.  3  but  for  Un  =  -100  m/s  and  L  =  100  m.  (Calculation  M) 
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GAMMA (SEC- I) 


M  =  1000.,  EL  =  20. M,  U0  =  -20. M/S 


Fig.  5.  As  in  Fig.  3  but  for  Un  =  -20  m/s  and  L  =  20  m.  (Calculation  S) 
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120  mesh  is  stretched  away  from  the  center  in  x,  where  the  grid 
size  is  10  m  in  both  directions.  A  small  initial  random 
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Fig.  7.  Plasma  configuration  at  10.67  sec  corresponding  to  the  initial 
conditions  in  Fig.  6. 
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perturbation  is  fifty  times  larger. 
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LAUREL,  MD  20810 

01CY  ATTN  DOC  LIBRARIAN 

01CY  ATTN  THOMAS  POTEMRA 

01CY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORPORATION 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
01CY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED 
STUDIES 

816  STATE  STREET 
(P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 


01CY 

ATTN 

DASIAC 

01CY 

ATTN 

WARREN  S.  KNAPP 

01CY 

ATTN 

WILLIAM  MCNAMARA 

01CY 

ATTN 

B.  GAMBILL 

LINKABIT  CORPORATION 

10453  ROSELLE 

SAN  DIEGO,  CA  92121 

01CY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 

01CY  ATTN  DEPT  60-12 

01CY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  8.  SPACE  CO .  ,  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 


01CY 

ATTN 

MARTIN  WALT/ 

DEPT  52-12 

01CY 

ATTN 

W.L.  IMHOF/ 

DEPT.  52-12 

01CY 

ATTN 

RICHARD  G.  JOHNSON/ 
DEPT.  52-12 

01 CY 

ATTN 

J.B.  CLADIS/ 

DEPT.  52-12 
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MARTIN  MARIETTA  CORPORATION 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

Old  ATTN  R.  HEFFNER 

MCDONNELL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  02647 
01CY  ATTN  N.  HARRIS 

Old  ATTN  J.  MOULE 

01CY  ATTN  GEORGE  MR02 

01CY  ATTN  V.  OLSON 

01CY  ATTN  R.V.  HALPRIN 

01CY  ATTN  TECHNICAL  LIBRARY 

SERVICES 


PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD 
LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD,  JR 
PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
UNIVERSITY  PARK,  PA  16802 
01  CY  ATTN  IONOSPHERIC 
RESEARCH  LAB 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  HA  01801 

01CY  ATTN  IRVING  L.  KOFSKY 


MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  03101 


Old 

ATTN 

P.  FISCHER 

Old 

ATTN 

W.F.  CREVIER 

Old 

ATTN 

STEVEN  L.  GUTSCHE 

Old 

ATTN 

R.  BOGUSCH 

Old 

ATTN 

R.  HENDRICK 

Old 

ATTN 

RALPH  KILB 

Old 

ATTN 

DAVE  SOVLE 

Old 

ATTN 

F.  FAJEN 

Old 

ATTN 

M.  SCHEIBE 

Old 

ATTN 

CONRAD  L.  LONGMIRE 

Old 

ATTN 

B.  WHITE 

01CY 

ATTN 

R.  STAGAT 

01CY 

ATTN 

D.  KNEPP 

Old 

ATTN 

C.  RINO 

MISSION  RESEARCH  CORPORATION 
1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NM  87106 

01CY  ATTN  R.  STELLINGVERF 

Old  ATTN  M.  ALME 

Old  ATTN  L.  WRIGHT 

MITRE  CORPORATION 
VESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

01CY  ATTN  W.  HALL 

01CY  ATTN  W.  FOSTER 


PHYSICAL  DYNAMICS,  INC. 


P.O.  BOX 

10367 

OAKLAND, 

CA  04610 

01CY 

ATTN 

A.  THOMSON 

RAD  ASSOCIATES 

P.O.  BOX 

9695 

MARINA  DEL  REY, 

CA  90291 

01CY 

ATTN 

FORREST  GILMORE 

01CY 

ATTN 

W.B.  WRIGHT,  JR 

Old 

ATTN 

V.J.  KARZAS 

Old 

ATTN 

H.  ORY 

01CY 

ATTN 

C.  MACDONALD 

Old 

ATTN 

BRIAN  LAMB 

01CY 

ATTN 

MORGAN  GROVER 

RAYTHEON  CORPORATION 
528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

Old  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  4 2ND  STREET 
NEW  YORK,  NY  10036 

01CY  ATTN  VINCF.  TRAPANI 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
10260  CAMPUS  POINT  DRIVE 


SAN  DIEGO 

,  CA  92 

.121-1522 

01CY 

ATTN 

L.M.  LINSON 

01CY 

ATTN 

D.A.  HAMLIN 

01CY 

ATTN 

E.  FRIEMAN 

01CY 

ATTN 

E.A.  STRAKER 

01CY 

ATTN 

C.A.  SMITH 
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SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1710  GOODRIDGE  DRIVE 
MCLEAN,  VA  22102 

01CY  ATTN  J.  COCKAYNE 

01CY  ATTN  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSVOOD  AVENUE 


MENLO  PARK,  CA 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

01CY 

ATTN 

94025 
J.  CASPER 
DONALD  NEILSON 
ALAN  BURNS 
G.  SMITH 
R.  TSUNODA 
D. A.  JOHNSON 
tf.G.  CHESNUT 
C.L.  RING 
VALTER  JAYE 
J.  VICKREY 
R.L.  LEADABRAND 
G.  CARPENTER 
G.  PRICE 
R.  LIVINGSTON 
V.  GONZALES 
D-  MCDANIEL 


TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  V.P.  BOQUIST 

TRW  DEFENSE  &  SPACE  SYS  GROUP 


ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 


Old 

ATTN 

R.K.  PLEBUCH 

Old 

ATTN 

S.  ALTSCHULER 

Old 

ATTN 

D.  DEE 

Old 

ATTN 

D.  STOCKVELL/ 

SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 

BURLINGTON,  MA 

01803 

01CY 

ATTN 

V.  REIDY 

01CY 

ATTN 

J.  CARPENTER 

01CY 

ATTN 

C.  HUMPHREY 

UNIVERSITY 

'  OF  PITTSBURGH 

PITTSBURGH 

!,  PA 

15213 

01CY 

ATTN 

N.  2ABUSKY 

